Phantom and non-phantom dark energy: The cosmological relevance of non-locally 

corrected gravity 
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In this paper we have investigated the cosmological dynamics of non-locally corrected gravity 
involving a function of the inverse d'Alembertian of the Ricci scalar, f(\J~ 1 R)). Casting the dy- 
namical equations into local form, we derive the fixed points of the dynamics and demonstrate 
the existence and stability of a one parameter family of dark energy solutions for a simple choice, 
/(□ _1 7?) ~ exp(aD _1 _R). The effective EoS parameter is given by, w c g = [a — l)/(3a — 1) and 
the stability of the solutions is guaranteed provided that 1/3 < a < 2/3. For 1/3 < a < 1/2 and 
1/2 < a < 2/3, the underlying system exhibits phantom and non-phantom behavior respectively; 
the de Sitter solution corresponds to a — 1/2. For a wide range of initial conditions, the system 
mimics dust like behavior before reaching the stable fixed point. The late time phantom phase is 
achieved without involving negative kinetic energy fields. A brief discussion on the entropy of de 
Sitter space in non-local model is included. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

At present, there are two major theoretical approaches 
to late time acceleration of universe which is supported 
by different data sets of observations. The standard lore 
is related to the modification of right hand side of the Ein- 
stein equations supplementing the stress energy tensor by 
a dark energy component Recently, serious attempts 
have been made to modify the geometry itself or the orig- 
inal Einstein-Hilbert action. A large number of papers 
are currently devoted to the investigations of f(R) gravi- 
ties (see review Q and references therein). These theories 
are motivated by phcnomcnological considerations. The 
problems faced by these theories can be circumvented 
in specific modelsQ. It is really interesting that these 
models do not reduce to cosmological constant A in the 
low curvature regime and thus can be distinguished from 
the latter. However, in any proposal on modification of 
gravity at large scales, it becomes mandatory to check 
whether the local gravity constraints are satisfied. The 
latter can be achieved for specific f(R) gravity models 
provided one invokes the chameleon scenario [3, 0| a la- 
Greek epicycle. 

It is also of crucial importance to explore the possibil- 
ity of obtaining late time acceleration from a fundamental 
theory. Thus the string curvature corrections to gravity 
and their cosmological relevance is a subject of current 
interest. Attempts have recently been made to derive 
current acceleration using the Gauss-Bonnet (GB) term 
and higher order curvature invariants coupled to a dy- 
namically evolving scalar field ||. The model with GB 
invariant exhibits a remarkable property that it does not 
disturb the scaling regime and can give rise to late time 
transition from matter regime to late time acceleration 
(see Ref.Qand references therein). Unfortunately, the 



coupling of GB term to scalar field gets large at late 
times; the nucleosynthesis also imposes stringent con- 
straints on these models. Inclusion of higher order cor- 
rections introduces further technical complications 0- 

Most of the proposals aimed to describe the late time 
cosmic evolution are faced with one or the other prob- 
lem. Recently, an interesting idea of using non-locally 
corrected Einstein theory was put forward in [1,0 . These 
corrections typically involve combinations of inverse of 
d'Alembertian of the Ricci scalar and might be induced 
by quantum loops and (or) stringy considerations. Be- 
ing non-local in character, these extra terms added to 
Einstein-Hilbert action can lead to late time accelera- 
tion as a time delayed effect avoiding the fine tuning 
problem Q. The non-local dynamics can be cast in a 
local form by introducing a number of auxiliary fields 0. 

In this paper, we consider the simplest form of non- 
local corrections to gravity and investigate the underlying 
cosmological dynamics in details. We explore the possi- 
bility of a matter like regime which can finally mimic 
(phantom) dark energy and in particular the de Sitter 
solution, at late times. The entropy of de Sitter space in 
non-local gravity is also briefly discussed. 



II. NON-LOCAL COSMOLOGY AND ITS LATE 
TIME ATTRACTORS 



In what follows we consider the local form of non- 
locally corrected gravity. Let us consider the following 
simple example of the non-local action Q 



S = 



(i) 



2 



where / is some function of d'Alembertian and denoted 
by □. The action would lead to non-local equations of 
motion which are difficult to investigate. The above ac- 
tion can be cast into local form by introducing two scalar 
fields and £@] 



S = 



d A x, 



= d xy^g 



+ /(<£)) + 



(2) 



The equivalence of actions |T]) and © can easily be 
checked. Indeed, by varying the action over £, we ob- 
tain 



U<j> = R or <j> = U~ l R, 



(3) 



which after substitution into ([2]) , gives the original action 
©. 

The local equations of motion can now be derived by 
varying with respect to the metric and the auxiliary fields 
<t> and £. Variation of ^ with respect to g^ u gives 

o = ^{iz(i + /(0)-O-fipW} 

-R„u (i + f{4>) -0 + 1 id^d v 4> + d^dvt) 
- (g^D - v M v„) Ui4>) -0 + ^T^ . (4) 

The variation with respect to <f> leads to 

= + (5) 

For obvious reason, we shall specialize to FRW back- 
ground with the metric 



ds 2 = -dt 2 



i(t) 2 (dx l ) 



i=l,2.3 



in which case the scalar fields an d £ are functions of 
time alone. In the FRW background, the time-time and 
the space-space components of Eq. @ yield the modified 
Hubble equation and equation for acceleration 



-m 2 (i + - o - 3H ( f 







+7^0 + K 2 p , 

= (2H + 3H 2 ) (1 + /(«/,)-£) + 1^6 



dt 



(7) 



(8) 



where p and p refer to the energy density and pressure 
of the background fluid and satisfy the usual continuity 
equation 



p + 3H(l + w)p = Q 



(9) 



with w being the equation of state parameter of the back- 
ground matter. The scalar field equations assume the 
following form 

= + 3Hj) + 6H + 12H 2 , (10) 

= £+3#£- (6H+ 12H 2 ^ f(4>) . (11) 

For the sake of simplicity, we shall use the exponential 
form of /(</>) 



(12) 



where fo and a are constants. We shall now cast the 
underlying dynamical equations in the autonomous form 
and look for dark energy solution as stable fixed points 
of the dynamics. We use the following notations (hence- 
forth, we use the unit k 2 = 1) 



M 1 ^ ( 



(13) 



to arrive at the autonomous form of evolution equations 

(14) 



dx H , 

dN = -12 --(., + 6) 



J| = -3y+12a-axy- j^(y ~ 6a) , (15) 



dz (n o o H \ 



and the constraint equation 



i 



G 



xy + ax — y 



(16) 



(17) 



Differentiating Eq. ([7]) with respect to time and using the 
continuity equation ([§]), we can express H / H 2 in the term 
of the autonomous variables as 



(6) = 



H 24a + Aax - a 2 x 2 - (4 + x)y - (1 + w)z 
H 2 ~ = 2[(l + f)(-6a + y) + f] 



(18) 



First we shall analyze the dynamics in the absence of 
the background fluid. In this case, the system (fl4]) — (jT5|) 
possess two critical points, (x c , y c ) — (0, 0) and (x c , y c ) — 
(—2/ a, 6a/ (2 — 3a)). The stability of the critical points 
can easily be analyzed. The eigenvalues of the stability 
matrix for the autonomous system corresponding to the 
first critical point (x c ,y c ) = (0,0) are given by, { — 1, — 1} 
making it a stable node. From Eq. lfTB)) , we easily find 
the effective EoS parameter at the critical point under 
consideration 



Wcff 



2_R_ 

3H 2 



(0,0) 



(19) 



For the fixed point (x c ,y c ) = (— 2/a, 6a/(2 — 3a)), using 
Eq.®, we find 



a- 1 
3a — 1 



(20) 
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FIG. 1: The evolution of w e a from the different initial con- 
ditions, dark black line corresponds to a — 0.5, x = —1.99, 
y = 5322.02, gray to a = 0.487, x = -1.99, y = 1012.06 and 
dashed to a = 0.507, x = -1.99,2/ = 1013.06. At the fixed 
point w e g= —1, —1.11 and —0.95 respectively 
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FIG. 2: The phase portrait of the system in absence of matter 
for a = 1/2. The phase space splits into two disjoint regions: 
all the trajectories starting from y > 3 converge to de Sitter 
where as they approach (0,0) for y < 3 and x > —4. 



In this case, the eigenvalues of the stability matrix are 
given by, (3a - 2)/(3a - 1), (3a - 2)/(3a - 1). Thus the 
fixed point is stable provided that l/3<a<2/3 leading 
to dark energy solutions, — oo < w e g < —1/3; the stable 
de Sitter solution is obtained for a = 1/2, see FigJTJ In 
FigUl we have plotted the phase portrait demonstrating 
the stability of both the critical points. Without the loss 
of generality, we have assumed both <fi and <f> to be nega- 
tive which corresponds to a particular type of boundary 
condition in Eq.©[||. Inclusion of matter in the system 
makes the analysis cumbersome in general. However, one 
can easily obtain the fixed point for H = 



12a 
4a- 3 ' 



8a(2a- 1) 
1 + w 



(21) 



In this case, the eigenvalues of the stability matrix are 

{-3,4a- 3,4a- 3- 3w} . (22) 

with a restriction that a = 3 (1 + w) /4. In this case, the 
fixed point is stable provided that w < 0. 



Despite the fact that it is extremely difficult to treat 
this system analytically, several exact solutions may be 
found. For instance, in case of the de Sitter solution, 
Eq. (fTU|) can be solved as 

3 -3H t , x 



-4H t 



(23) 



where (j)Q and <p\ are the constants of integration. For 
simplicity, we assume that 0o = <t>\ — 0. Using Eq. lfTTj) . 
we get the expression for £, 



t = 3 /o c -4qg t _ & c -3H t _ t 

6 4a -3 3H 



(24) 



Here £o and £i are constants. In the absence of matter, 
p = 0, the de Sitter space corresponds to a = 1/2 and 
£i = —1. When p 7^ 0, there exists a de Sitter solution 
provided we choose 



K 2 p 



3iJ2 (l + 3w) 



6 = -1 . (25) 



In the presence of matter with w ^ 0, one may find a 
de Sitter solution even for a more complicated choice of 

m = foe^ 2 + Ae 3 ^ +1 )*/ 4 . (26) 
The solution is given by 

<j> = -4H t , f = 1 + 3 f e' 2Hot + A e -3(^+i)ffo t j 
3(3w + l)tf 2 /i 



-3{l+w)H t 



(27) 



Let us now briefly point out some generalities of our 
formulation. The above discussion of non-local cosmol- 
ogy is based upon the simplest version of non-locality. 
Indeed, we could start from a general non-local action: 

S = J [F (R, OR, D 2 R, • • • , D m R, D^R, 

D- 2 R, • • • , n-»R) + £ ma tter] ■ (28) 

with m and n being the positive integers. By introducing 
scalar fields A, B, Xk, Vk, (k = 1,2, • • • ,m), and 0/, 
(I = 1, 2, • • ■ , n), we may rewrite the action (|2"5|) in the 
following form: 



S = 



d 4 x. 



-g [BR — BA + F (A,r]i,r]2 7 - ■ ■ ,rj m , 

m 



'Ik- 1 



k=2 



1=1 

^^matter] ■ 



fe=l 



1=1 



(29) 



The variations over A, B \k, £z leads to the following 
equations 

Q = R-A = r )1 -nA = r) k - = D<pi — A 
= Ufa - 4,1-1 {k = 2, • • • , m, I = 2, • • • ,n) . (30) 
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which gives % = O k R and □ 0j = i? and we find that 
the actions (j2"9")l and (j2"5j) are equivalent. 

Let us now consider a scale transformation given by 



1 



(31) 



Then the action (|2"9")l is transformed into the action in 
the Einstein frame: 



S= d 4 



1 / m n \ 

+ 2£ + £ ^Xkd^Vk-i + E W>i 

V fc=2 (=1 / 



matter 



fc=l 



;=2 



(32) 



Here £^ attor can be obtained by scale transforming of 



the metric tensor in £ ma tter by (|31[) . The above action 
can be regarded as a non-linear a model with potential 
coupled to gravity. Action (Q]) can be thought as a par- 
ticular case of the general non-local action. In this case, 
however, the continuity equation ^ is no longer valid; 
additional terms on the right hand side of this equation 
are generated and are responsible for inducing interaction 
of matter with background [4|. This aspect is crucial for 
the discussion of local gravity constraints on the models 
based on the modified theories of gravity. 



III. THE ENTROPY OF DE SITTER SPACE 

In this Section, we shall investigate the de Sitter space 
entropy properties in non-local modified gravity with the 
cosmological constant and matter. 

The starting point is the trace of equation of motion of 
the non-local model, which reads, for a constant curva- 
ture solution and in the presence of a cosmological con- 
stant A and matter, as 

C/W + 1-0^ = 4A-k 2 T 

+6f((t>)R + 3f"(ct>)d p ct>d>>ct> +d p tdr<t>, (33) 

with T the stress tensor trace. Let us suppose that 
spherically-symmetric, static, constant curvature solu- 
tion exists 



ds 2 = -A(r)dt 2 s + Air^dr 2 + r 2 dY? 



(34) 



There is a horizon if A(th) = with A'(rjj) 7^ and 
rjj > 0. Then there is also an entropy associated with 
this horizon, entropy that may be evaluated using Wald's 
method. For the local action, one has 



S. 



BH 



Ah 
AG 



[1-C + /(0)L 



(35) 



with all the fields being evaluated on the horizon and 
Ajj = Airr H . Thus, in general, there is a violation of the 
"Area Law" . 

Within the above black hole like metric, Eq. D<j) = 
R and Eq. for the field £ read for time dependent and 
spherically symmetric fields as 



d 2 



1 



.4 



d r (r 2 Ad r 



-^ + ^dr(r 2 Ad r £) 



R, 



-f'MR. (36) 



One may look for the solution of the first equation in the 
form (j> = — 4-ffo^s — 2 In A(r). As a result, it follows that 
de Sitter space Ads(r) = 1 — Rr 2 / 12 is a solution. Then, 
it is easy to show that the solution for £ reads 



&r,t e )=B- 



3/o 



(4a - 3) 



(1 - H 2 r 2 )~ 2a e 



-AaH t, 



(37) 



As a result, the entropy factor becomes 
(4a-6)/ , 



!-£+/(<£) = 1-B- 



(4a - 3) 



(38) 

As a consequence, for a > 0, the field £ and the entropy 
are, in general, divergent on the horizon, while, as soon 
as a < 0, the £ field and the de Sitter entropy are finite 
and independent on the time t s on the horizon. However, 
there exists a positive value, namely a = 3/2, which ren- 
ders the de Sitter entropy finite. We will see the physical 
significance of these two choices later. 

Furthermore, in order to satisfy Einstein equations, the 
choice made for /(</>) has to satisfy the trace constraint 
(J33J) , where the stress tensor of the matter and cosmolog- 
ical constant contribution are present. A direct computa- 
tion shows that the trace constraint is satisfied as soon as 
1— B = A/3H 2 , and the stress tensor trace is not zero and 



is given by T = T A ds {r)- 2a e 



2a a -4H at. 



, where Tq depends 



on a. Again, this trace diverges on the horizon as soon as 
a > 0. Note that in the case a = |, the entropy of the de 
Sitter solution reads Sds = (Ah /4G) (A/3H 2 ). For van- 
ishing cosmological constant we have found a de Sitter 
solution with finite entropy which is vanishing! However, 
recall the Wald method gives the entropy modulo a con- 
stant. 

In the FRW space-time, the solutions we have obtained 
lead to 



A 



K 2 (l-W)p 3Ha (l + w )t 



3H 2 3H 2 w(l + 3w) ' 
while equation of motion for £ gives 

4a - 6 



1 - e + m = i-B+f c 

A comparison with ui^l gives 
A 



-4aH t 



4a 



, (39) 



(40) 



1 - B = 



3H 2 



4a = 3(1 + ib), f 



k 2 p 



3H 2 (1 



3w) ' 
(41) 



5 



while if w = 1, we have fo undetermined. As a check, we 
may pass to the static de Sitter patch by transforming 
the time coordinate t — t s + (\/2Hq) ln(l — H^r 2 ), and 
note that one gets equation ([55]) again. 

Thus, the existence conditions found for the de Sit- 
ter entropy require a non vanishing positive cosmological 
constant and matter and, in general, this matter, in or- 
der to have finite entropy, should be "phantom" matter, 
since the condition a < is equivalent to 1 + w < 0. 
However, as already noted in the static patch, if we have 
matter such as w = 1, namely stiff matter p = p, im- 
portant in the early universe, the time dependence of 
the entropy drops out, in agreement with the result ob- 
tained in the static gauge. One can also see that entropy 
divergences in the static gauge correspond to the (non- 
physical) time-dependence of entropy in the cosmological 
gauge and thus this situation seems problematic. This 
consideration indicates that perhaps even the definition 
of entropy for non-local gravity should be reconsidered. 
It is also interesting that AdS black hole solution for the 
theory under consideration does not exist, while de Sitter 
solution exists. 



IV. CONCLUSIONS 

In this paper we have investigated the cosmological 
relevance of non-local corrections to gravity of the type 
~ exp(aD _1 i?) which can be motivated by 
non-perturbative quantum effects. Casting the equations 
in local form we analyzed the underlying dynamics and 
explored the possibility of finding stable dark energy so- 
lutions. We have shown that the system might mimic 
the matter dominated regime for a long time and ul- 
timately switches over to dark energy solutions at late 
times thereby alleviating the fine tuning problem. The 
class of the dark energy solutions is parameterized by a 
which ranges as 1/3 < a < 2/3. The corresponding ef- 
fective EoS parameter varies as, — oo < w c ff < —1/3. In 
case of a > 1/2, we have non-phantom dark energy; the 
de Sitter solution is obtained for a = 1/2. It is remark- 
able that there exists a range of the parameter, namely, 
1/3 < a < 1/2 for which the model leads to phantom 



dark energy solutions which are attractors of the system. 
It should be emphasized that the observed acceleration 
can be achieved for natural values of a ~ 1. Secondly, 
the model under consideration does not involve negative 
kinetic energy fields. Unlike standard phantom theories 
giving rise to transient phantom phase, the non-locally 
corrected gravity leads to a class of phantom energy solu- 
tions which are the late time attractors of the dynamics. 

It is also important to note that even if more pre- 
cise observational data defines the EoS parameter to be 
slightly different from —1, there exists a possibility of re- 
alizing such scenario a in non-local gravity as the effective 
quintessence or phantom cosmology. 

Interestingly, the existence of de Sitter entropy in the 
non-local model requires a non vanishing positive cosmo- 
logical constant and phantom matter. 

Last but not the least, we should remember that the 
low energy modifications of gravity crucially differ from 
the corresponding local versions. In general, any large 
scale modification of gravity faces two major challenges: 
The first is related to the presence of either the ghost or 
the tachyon modes in these models; secondly, the solar 
physics imposes stringent constraints on these models. 
Thus it would be interesting to carry out these investi- 
gations for non-local theories, as discussed in this paper; 
we defer this analysis to our future work. It would also 
be interesting to investigate the complicated versions of 
non-local gravity mentioned at the end of section II. 
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